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1.  INTRODUCTION 


Radioactive  argon  gas  (A-41)  is  produced  when  naturally  occurring 
A-40  captures  a neutron.  Since  air  contains  approximately  1 percent 
A-40,  neutron  irradiation  of  air  produces  a measurable  quantity  of  A-41, 
which  decays  with  a half-life  of  110  min. 

At  the  Diamond  Ordnance  Radiation  Facility  (DORF) , a research 
reactor  is  operated  adjacent  to  a large  air-filled  volume  called  the 
exposure  room.  The  irradiated  air,  when  exhausted  from  the  room  through 
a gas  stack,  contains  a measurable  quantity  of  A-41.  To  be  certain  that 
this  radioactivity  does  not  represent  a radiation  hazard  to  the 
surrounding  environment,  the  release  of  A-41  must  be  measured  and 
reported  annually  as  required  by  Army  Regulations  385-80. 

To  assist  in  the  analysis  of  the  amount  of  argon  released,  a 
computer  program  has  been  written  to  take  the  data  recorded  on  the 
gas-stack  monitor  and  calculate  the  conversion  constant  between  counts 
per  minute  on  the  monitor  and  the  number  of  curies  of  A-41  released  to 
the  environment.  The  procedure  for  obtaining  the  experimental  input 
data  is  given  in  a memorandum  by  T.  Chichester.1  This  report  documents 
the  procedure  used  to  process  and  record  the  data  for  the 
gas-stack  monitor  calibration. 


2 . BACKGROUND 

The  DORF  reactor  is  a TRIGA-type  reactor  capable  of  operating  in  a 
pulsed  or  steady-state  mode.  The  reactor  is  movable  and  can  operate 
either  completely  shielded  by  water  or  immediately  adjacent  to  an 
exposure  room.  During  operation  next  to  the  exposure  room,  air  (which 
is  flowing  through  the  room)  is  irrradiated  and  vented  to  the  outside 
through  the  gas  stack  (fig.  1,  2,  3).  Several  radioactive  effluents  are 
produced;  however,  A-41  is  the  only  important  one,  because  the  others 
have  sufficiently  short  half-lives  that  they  decay  before  leaving  the 
stack. 

Located  near  the  lower  end  of  the  air  stack  is  the  gas-stack 
monitor.  This  system  (fig.  4)  continuously  samples  the  air  in  the 
stack.  The  sampled  air  is  pumped  into  a chamber  containing  a Geiger 
tube.  A logarithmic  circuit  accepts  this  output  voltage  and  converts  it 


Procedure  for  Calibrating  the  Gas-Stack  Monitor , Health  Physics 
Section , Walter  Reed  Army  Medical  Center . 
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Figure  1.  Diamond  Ordnance  Radiation  Facility,  main  building 
with  a view  of  the  gas  stack. 


CARRIAGE 

assembly 


Figure  2.  Vertical  section  of  the  DORF  reactor. 
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Figure  3.  Perspective  view  of  DORF  reactor. 


to  a logarithm  output  that  drives  a count  rate  meter  and  a strip-chart 
recorder.  It  is  necessary  to  relate  the  trace  of  the  recorder,  in  units 
of  counts  per  minute,  to  the  number  of  curies  per  cubic  centimeter  in 
the  air  sampling  chamber.  When  this  relationship  is  known,  the  number 
of  curies  per  minute  released  is  found  by  multiplication  with  the  stack 
flow  rate. 

The  process  for  finding  the  conversion  constant  between  counts  per 
minute  and  curies  per  cubic  centimeter  begins  by  irradiation  of  a 210-cc 
flask  of  99.97  percent  pure  argon.  This  irradiation  produces  a high 
concentration  of  A-41.  A second,  unirradiated  flask  of  A-40  is  placed 
in  the  gas-stack  monitor  system  in  series  with  the  irradiated  flask  by 
means  of  tubing  (fig.  5) . The  contents  of  both  flasks  are  pumped 
through  the  gas-stack  monitor  by  a varistaltic  pump  to  allow 
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r±gure  4.  Gas-stack  monitor. 


equalization  of  the  A-41  concentration  throughout  the  system.  When 
equilibrium  is  attained,  the  irradiated  flask  is  sent  to  the  Health 
Physics  Office,  Walter  Reed  Army  Medical  Center,  for  radio-analysis. 
The  other  flask  is  left  connected  to  the  gas-stack  monitor. 
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IRRADfATED  FLASK 


UN  IRRADIATED  FLASK 


Figure  5.  Diagram  of  the  argon-41  flow  system  during  calibration. 


The  sample  is  analyzed  at  the  Health  Physics  Office.  A report  sent 
to  DORF  gives  the  concentration  of  A-41  in  the  flask  in  microcuries  per 
cubic  centimeter  at  a given  time.  After  this  concentration  is  corrected 
for  decay  to  a predetermined  time,  it  is  possible  to  correlate  the  count 
rate  on  the  gas- stack  monitor  with  the  concentration  of  the  A-41  in  the t 
sample  chamber. 

3.  DATA  PROCESSING  PROCEDURE 

3. 1 Handling  the  Raw  Data 

When  the  unirradiated  flask  is  sent  to  the  Health  Physics 
Office,  the  gas-stack  monitor,  still  connected  to  the  irradiated  flask, 
is  left  running  for  about  18  hr.  In  this  time,  the  logarithmic  trace  of 
the  recorder  comes  on  scale  at  100,000  cpm  and  goes  through  a number  of 
half-lives  before  it  blends  with  the  background  activity  (fig.  6) . A 
straight  line  is  drawn  from  the  point  at  which  the  trace  comes  on  scale 
to  the  low  count  rate  end  where  the  trace  starts  to  become  nonlinear 
(because  of  the  effects  of  background) . If  the  line  has  the  correct 
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Figure  6.  Trace  of  the  gas-stack  monitor  data. 
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slope  (i.e.,  exhibits  the  110-min  decay  of  A-41),  after  the  data  have 
been  fitted  to  a straight  line  by  the  method  of  least  squares,  there  is 
no  indication  of  a leak  in  the  system,  and  the  calibration  is  assumed  to 
be  valid. 


If  the  slope  is  not  correct,  the  cause  must  be  found,  and  the 
procedure  must  be  repeated.  Possible  causes  are  leaks  in  the  system  and 
inadequate  background  compensation.  The  time  at  which  the  trace  comes 
on  scale  is  defined  as  time  zero.  The  sample  from  the  Health  Physics 
Office  is  decay  corrected  to  time  zero.  Then  starting  at  time  zero,  the 
count  rate  is  taken  at  110-min  intervals  (i.e.,  every  half-life),  and 
the  ratios  of  A-41  concentration  to  the  count  rate  are  calculated. 
These  data  are  taken  until  nonlinearity  becomes  a problem.  Both  the 
average  value  of  the  ratios  and  their  RMS  deviations  are  determined  for 
the  first  four  (i.e.,  most  statistically  significant)  data  points. 

3. 2 Details  of  Data  Processing 

The  computer  program  ARGON  that  performs  the  analysis  is 
currently  stored  in  the  DORF  library  on  time-share  computer  at  the 
Applied  Physics  Laboratories  in  Columbia,  MD.  The  program  (app  A)  is 
written  in  the  Conversational  Programming  System  (CPS) . 

To  execute  the  program,  the  user  must  first  call  the 
remote- terminal  port  of  the  computer  and  log  on  with  a valid  account 
number.  When  the  computer  answers  with  "READY,"  the  directive  "CPS"  is 
given,  telling  the  computer  what  working  language  is  being  used.  Then, 
the  user  can  begin  processing  data  through  the  teletype  in  accordance 
with  the  following  paragraphs.  An  example  is  given  in  appendix  B. 

3.2.1  Input 

To  begin,  the  program  is  transferred  from  the  user 
library  with  the  LOAD  (ARGON)  command.  A question  mark  follows,  and  the 
user  types  "XEQ."  Execution  of  ARGON  begins  at  this  time. 

The  first  input,  to  be  entered  after  "DATE(l)?"  is  typed, 
is  the  date  on  which  the  calibration  was  made,  entered  as  day,  month, 
and  year.  Next  (after  "NUM?"  is  printed) , the  number  of  half-lives 
measured  is  entered.  In  the  example  shown  in  appendix  A,  eight 
half-lives  were  measured.  The  concentration  of  A-41  in  the  assayed 
flask  is  the  next  input  parameter.  The  concentration  in  microcuries  per 
cubic  centimeter,  decay  corrected  to  time  zero,  is  entered  following 
"MICRO?."  The  remaining  values  to  be  entered  are  the  count  rate  from 
zero  time,  "COUNT  (0),"  to  the  last  half-life  value,  which  is  "COUNT 
(8)"  in  the  example.  The  count  rate  at  zero  time  is  always  100,000  cpm 
with  the  present  system,  and  the  remaining  count  rates  are  found  by 
measurement  of  110-min  intervals  from  zero  time  (fig.  6) . 
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3.2.2  Output 


The  analysis  (app  C)  follows  the  input  of  count  rates. 
All  input  data  are  listed , as  well  as  the  following  information: 

a.  A determination  of  the  half-life  and  count  rate  at  zero  time  by 
the  method  of  least  squares 

b.  The  error  in  the  half-life 

c.  The  concentration  of  A-41  at  each  half-life 

d.  The  calculated  ratio  of  microcuries  per  cubic  centimeter  per 
counts  per  minute  at  every  half-life 

e.  The  average  ratio  of  the  first  four  values  given  in  d 

f.  The  RMS  variation  of  the  first  four  ratios. 

The  output  shown  in  appendix  C is  from  the  calibration 
run  of  27  August  1974.  Results  show  that  there  were  no  leaks  in  the 
system;  however,  results  show  also  inadequate  background  compensation. 
This  is  indicated  by  an  increase  in  the  calculated  half-life  by 
3.5  percent  after  five  half-lives.  The  low  RMS  error  shown  indicates 

-the  value  of  the  better  statistics  of  the  first  four  points. 
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APPENDIX  A. --ARGON  PROGRAM 


10. 

/* 

ARGON 

♦ /; 

1 5. 

/*  AN  ARGON-41  ANALYSIS  PROGRAM  */; 

E0  . 

DECLARE  DATE  C 3),COUNTC0; 10); 

Zb. 

DECLARE  DDC0: 1 0), DRCu: 1 0) ! 

30. 

SUMA 

,jSumxy*sumae=G; 

3b. 

SUMDEL, SUMRAT=0! 

40  • 

PUT 

LISTC  'ENTER  THE 

DAY, MONTH, YEAR  ' ) ; 

Ab. 

GET 

LIST(DMTE); 

50. 

PUT 

LISTC  1 ENTER  THE 

NUMBER  OF  Half  LIVES  MEASURED')! 

55. 

GET 

LiSTCNUM); 

60  • 

PUT 

LISTC  • ENTER  THE 

ASSAYED  CONCENTRATION  OF  ARGON-41  I 

N* 

) i 

65. 

PUT 

LISTC  'MICRU-CUR IES/CC  CORRECTED  I OR  THE  110  MIN.  HA 

>; 

70. 

PUT 

LISTC 'LIFE  TO  TIME  = 0 OF  THE  SE.KIES  OF  COUNTS')! 

75. 

GET 

LISTCMICRO)! 

80. 

PUT 

LISTC  'ENTER  THE 

INITIAL  COUNTS  AT  ZERO  HALF  LIVES') 

J 

85. 

GET 

LISTCCOUNTC0 ) )i 

90. 

PUT 

LISTC  'FIN  1 ALL Y , 

ENTER  THE  COUNTS  AT  EACH  FOLLOWING 

AL.F  LIFE*); 

95. 

T = 1 1 

0; 

100.  T 1 ME  = 0! 

105.  N=NUM+ 1 ! 

110.  SUMY  =LOGCCOUNTC0) ) ! 

115.  SUMY2=SUMY  *SUMY  J 

120.  DO  1=1  TO  NUM! 

125.  GET  LIST (COUNT Cl))! 

130.  T 1 ME  = T I ME  + T J 

135.  Y=LOGCCOUNTC I ) ) ! 


1 40  . 
1 45. 
1 50  . 
1 55- 
1 60  . 
1 65. 
1 70  • 
175. 
180. 
185. 
1 90  . 
1 95. 
200  . 


SUMX=SUMX+TlME! 

S UM Y =S  UMY  + Y J 
SUMXY=SUMXY*TiME*Yi 
SUMX2=SUMX2*TiME*TlME! 
SUMY2=SUMY2+Y*Y! 

END'  ! 

DEN=SUMX**2 -N*S  UMa2 ! 

A 1 =CSUMX*SUMY-N*SUMXY)/DEN! 

TEMP  = - 1 * CSUMX2*SUMY-SUMX*SUMXY )/DEN! 
AO  = EXPC  TEMP ) J 
DO  1=0  TO  NUM! 

D=MICR0/2**I ! 

KAT 1 0 = D/COUNT ( 1 )! 


? 

205.  SUMRAT=SUMRAT*RATIO! 

210.  DDC I ) =D* 1 0000; 

215.  DR  C I ) =RAT I 0 * • 1 E09! 

220.  END  ! 

225.  X8AR=SUMRaT* . 1 E09/N1 

230.  PUT  EDITC ' ')CSKIPC8),FC0))! 

235.  PUT  L I ST ( ’ LEhST  SQUARES  ANALYSES  OP  ARGON-41  CALI BRATI 0 

N*  ) ; 

240.  PUT  L 1 S T ( * 
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APPENDIX  A 


- ' )i 


245. 

2 50. 

255. 

260. 

C PER  C PM • ) J 
265. 

270. 

275. 

260. 

26  5 • 

2*0. 

295. 

300. 

B'  ) J 
305. 

310. 

31  5. 

320. 

32  5. 

330. 


PUT  LISTC  # # ) J 

PUT  LISTC  #DATE=  DATEC 1 ) # '/ ' > DA TEC  2 ),  '/  * » DATEC  3 ) > J 
PUT  LISTC  ' ' )i 

PUT  L I ST ( 'HP  VALUE  OF  CONCENTRAT I ON  = % MICRO,  # M IGROC I . /C 
PUT  LISTC • * ) J 

PUT  LISTC 'INPUT  DATA  ARE  AS  FOLLOWS:'); 

PUT  LISTC 'NO.  TIMECMIN)  CO UNTS/M IN . ' ) ; 

DO  1=0  TO  num; 

PUT  I MAGEC I * T* I > COUNTC I))ClME)i 
END  i 

PUT  LISTC  ' ' )i 

PUT  LISTC  'SLOPE  AND  INTERCEPT  OF  ARGON-41  DECAY  CURVE  A 

PUT  LISTC  ' ' ); 

PUT  1 MAGEC -LOGC  2 ) /A 1 )CIMB)J 

PUT  LISTC ' ' )J 

PUT  IMAGECAO) C IMC); 

PUT  LISTC ' ' )i 

ER=£BSC  C LOGC 2 ) /A1 +T)/T)*100J 


335. 

PUT 

I MAGEC  ER) 

3 40  . 

PUT 

LISTC  * * ) 

345. 

PUT 

LISTC  * * ) 

3 50. 

PUT 

LISTC  'CAL 

OWS : ' ) ; 

355. 

PUT 

LISTC  '*** 

*#.***  • ) j 

3 60. 

PUT 

LISTC  ' ' ) 

365. 

PUT 

LISTC'  ') 

370. 

PUT 

LISTC  * 

nos  • ); 

375. 

PUT 

LISTC  ' 

10-8)  * )) 

380. 

PUT 

LISTC  ' ' ) 

PUT  LISTC  ' it***************************************  **** 


OBSERVED 

COUNTS 


CALCULATED  ACTIVITIES  RA 
(X  10-4)  (X 


385. 

390. 

395. 

400. 

405. 

410. 

415. 

420. 

42  5. 

430  . 

435. 

440  . 

445. 

450. 

455. 

4 60.  I MAS 


DO  1=0  TO  NUM) 

PUT  I MAGEC COUNTC I >» DDC lit DRC 1 ) ) C IMA) J 
END  ) 

XBAR  = CDRC0)«-DRC  1 )+DRC2)+DRC3)  )/4) 

DO  1=0  TO  31 

SUMDEL  = SUMDEL«-CABAR-DRC  I ) >**2) 

END  ) 

RMS  = SQRT  C SUMDEL/3 ) /XBAR) 

PUT  LISTC  * * )) 

PUT  LISTC*  FOR  THE  MOST  SIGNIFICANT  DATA  POINTS*)) 
PUT  LISTC  * * ) ) 

PUT  IMAGECXBAR) C IMF)) 

PUT  LISTC  • * )) 

PUT  IMAGECRMS*100) C IMG)) 

PUT  EDITC  * * ) CSKIPC8)«FC0) )) 

IMAGE) 


465.  IMfct:  IMAGE) 
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APPENDIX  A 


CALCULATED  HALE  L I FE=  • - MIN 

470.  IMC:  IMAGE! 

INTERCEPT COUNTS 


475. 

IMD: 

IMAGE! 

ERROR  IN 

hale 

L I FE  = . - - 

- PERCENT 

480  * 

IME: 

IMAGE! 

485. 

IME: 

IMAGE! 

AVERAGE 

RATIO 

= MICROCURIES/CC/CPM 

490. 

IMG: 

IMAGE! 

R MS  VAR . 

IN  RAT  10  = • 

PERCENT 

495. 

STOP  I 
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APPENDIX  B. — ARGON  USER'S  GUIDE 


AEU 

ENTER  THE  DAY  t MONTH*  Y EAR 
DATEC 1 ) 

►27  « 74 

ENTER  THE  NUMBER  OF  HAlF  LIVES  MEASURED 
NUM 
► 8 

ENTER  THE  ASSAYED  CONC ENTRaT 1 ON  OF  ARGON-41  IN 
M1CR0-CUR1ES/CC  CORRECTED  FjR  THE  110  MIN.  HALF 
LIFE  TO  T1ME=0  OF  THE  SERIES  OF  COUNTS 
MICRO 

► .00t:  6b 

ENTER  THE  INITIAL  COUNTS  AT  ZERO  HALF  LIVES 
COUNT  C 0 ) 

► 1 00000 

^ INIAlLY * ENTER  THE  COUNTS  AT  EACH  FOLLOWING  HALF  LIFE 
COUNTC 1 ) 

► S0000 
COUNT (2  > 

►26000 
COUNT (3 ) 

► 12500 
COUNTC 4) 

► 6500 
COUNTC  5) 

►0500 
COUNTC  6) 

► 1 900 
COUNTC 7) 

► 900 

COUNTC tt) 

►450 
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APPENDIX  C. — ARGON  OUTPUT 


LEAST  SQUARES  ANALYSES  OF  ARGON-41  CALIBRATION 


DATE  = 27/8/74 

HP  VALUE  OF  CONCENTRATIONS  .00268  MICROCI./CC  PER  CPM 


1 NPUT 

DATA  ARE  AS 

FOLLOWS: 

NO. 

T 1 ME (MIN) 

COUNTS/MI 

0 

0 

1 00000 

1 

1 1 0 

50000 

2 

220 

2 6000 

3 

330 

12S00 

4 

440 

6500 

5 

5 30 

3 500 

6 

660 

1 900 

7 

770 

900 

8 

880 

450 

SLOPE  AND  INTERCEPT  OF  ARGON-41  DECAY  CURVE  ARE 

CALCULATED  HALF  L I FE=  113.9  MIN 

INTERCEPTS  9802  S COUNTS 

ERROR  IN  HALF  L I FE=  3.526  PERCENT 


CALIBRATION  CONSTANTS  FOR  THIS  RUN  ARE  AS  FOLLOWS  : 

************************************************** 


OBSERVED 

CALCULATED  ACTIVITIES 

RATIOS 

COUNTS 

(X  10-4) 

<X  10-8) 

1 00000 

26.8000 

2. 6800 

S0000 

13.4000 

2 . 6800 

26000 

6.7000 

2.5769 

12  500 

3.3500 

2 . 6800 

6500 

1 . 6750 

2.5769 

3S00 

0.8375 

2 . 3929 

1 900 

0.4187 

2.2039 

900 

0.2094 

2.3264 

450 

0. 1047 

2.3264 

FOR  THE  MOST  SIGNIFICANT  DATA  POINTS 
AVERAGE  RATIO=  2.654  MICROCURIES/CC/CPm 
RMS  VAR.  IN  RAT10=  1.942  PERCENT 
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